Engineering vortex rings and systems for controlled studies of vortex 
interactions in Bose-Einstein condensates 
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We study controlled methods of preparing vortex configurations in atomic Bose-Einstein con- 
densates and their use in the studies of fundamental vortex scattering, reconnection processes and 
superfluid sound emission. We explore techniques of imprinting closed vortex rings by means of co- 
herently driving internal atomic transitions with electromagnetic fields which exhibit singular phase 
profiles. In particular, we show that a vortex ring can be prepared in a controlled way by two 
focused co-propagating Gaussian laser beams. More complex vortex systems may also be imprinted 
by directly superposing simpler field configurations or by programming their phase profiles on op- 
tical holograms. This provides the controlled method of studying vortex reconnections in atomic 
superfluids. We analyze specific examples of two merging vortex rings in a trapped two-component 
87 Rb condensate. We calculate the radiated sound energy in the vortex ring reconnection process 
and show that the vortex relaxation and the re-distribution of sound energy can be controlled by 
the imprinting process. The energy is first concentrated towards the trap center and later emitted 
outwards as sound and transformed to surface excitations. As another such technique, we study cre- 
ating pairs of 2D point vortices in Bose-Einstein condensates using a 'light roadblock' in ultra-slow 
light propagation. We show how this can be used to study vortex collisions in compressible super- 
fluids and how these collisions result in energy dissipation via phonons and, sometimes, annihilation 
of vortex pairs. 

PACS numbers: 03.75.Kk,03.75.Lm,03.75.Mn,67.40.Vs 



I. INTRODUCTION 

One of the advantages of atomic Bose-Einstein conden- 
sates (BECs), as compared to more traditional quantum 
fluids, is the dramatic flexibility of experimental prepa- 
ration. Both the internal and the center-of-mass states 
of ultra-cold atoms can be manipulated and controlled 
in diverse ways using electromagnetic (em) fields. In this 
paper we explore the possibilities of using the present day 
atomic physics technology as a sophisticated state engi- 
neering tool, in order to construct highly non-trivial su- 
perfluid states in atomic many-particle systems. In par- 
ticular, such a 'field-theory engineering' may be useful in 
preparing topological defects, studying their interactions, 
and the decay of superfluid turbulence Q . We consider 
em field configurations which exhibit phase singularities 
in the field amplitude. When such a field is used to drive 
internal atomic transitions, the non-uniform phase profile 
of the em field amplitude may be imprinted on the mat- 
ter wave, resulting in a topologically non-trivial phase 
profile in the atomic condensate. We show how a reason- 
ably simple em field configuration of two phase-coherent 
Gaussian laser beams could be used to imprint a vor- 
tex ring in atomic BECs by simultaneously controlling 
the position, the orientation, and the radius of the ring. 
More complex vortex structures, such as knotted vortex 
lines, may be imprinted using optical holograms. 

Moreover, we study in detail the application of this 
technique to the imprinting of a pair of vortex rings 
and study their reconnection dynamics in a trapped two- 



component Rb condensate, using experimentally fea- 
sible parameters. Such reconnections have escaped ex- 
perimental observation in BECs to date. The two rings 
merge together and the resulting turbulent dynamics gen- 
erates an effective local dissipation in the gas, allowing 
the vortex configuration to relax to lower energy without 
an explicit damping term in the dynamics. This mech- 
anism is also analogous to the dynamical formation of 
a vortex lattice in a rotating zero-temperature BEC |2|. 
We qualitatively analyze the re-distribution of sound en- 
ergy in the trapped cloud due to the vortex reconnection 
processes. Here the excitation of Kelvin waves and the 
radiated sound energy can be controlled by changing the 
parameters of the em fields. Although vortex rings, as 
localized singular defects which are detached from su- 
perfluid boundaries, are interesting in their own right, 
they can also form a building block of energetically sta- 
ble localized particle-like solitons (3D Skyrmions) in a 
two-component 87 Rb BEC [HOHH, providing a link be- 
tween ultra-cold atom physics, elementary particles, and 
cosmology. 

As another method of engineering systems for study- 
ing fundamental vortex interactions, we consider the gen- 
eration of a 'light roadblock' [|| by abruptly distort- 
ing the ultra-slow light propagation 0> S 13 m a BEC. 
The strong nonlinear coherent light-matter coupling is 
here based on the method of electromagnetically-induced 
transparency (EIT) [lfj, which allows propagation of 
light pulses through the BEC due to the lack of absorp- 
tion. If the probe and the coupling beams are prop- 
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agating in orthogonal directions and the coupling field 
is varied quickly to zero over a short distance near the 
BEC center, the probe field is abruptly compressed and 
stopped inside the atom cloud. This results in large 
transfer of atom population between the internal states 
and the emergence of a very narrow density defect inside 
the BEC, which subsequently generates solitons via quan- 
tum shock waves p| and then vortex pairs via the snake 
instability [TlllllTl^|. In Ref. [3 it was noted that this 
can be used to create a 'gas' of multiple vortex particles, 
of both circulations, which are out of equilibrium and 
subsequently interact. Here we show that in a 2D geom- 
etry the generation of the defects can be controlled to 
yield desired configurations of point vortex pairs whose 
collisions and sound emission can be investigated in de- 
tail. In particular, we focus on a regime in which the size 
of the vortex pairs is comparable to the healing length in 
the superfluid. In such a case, collisions of vortex pairs 
are more complicated than they are in the imcompress- 
ible limit 0] , as the collisions can cause energy to be dis- 
sipated in the form of phonons. This dissipation alters 
the collision dynamics and can lead to annihilations of 
vortex pairs, another interesting phenomenon which has 
hitherto been unobserved experimentally. We present nu- 
merical study of this and show how the light roadblock 
may allow its observation. 

The quantized vorticity in atomic superfluids has in- 
spired considerable experimental and theoretical activity 
in recent years pl| . There have been several theoretical 
proposals to imprint vortex line singularities on atomic 
BECs by means of transferring angu lar momentum on 
atoms from em fields [13, [IM [la H3, HH H3, H3, 13, 
and some of the techniques have alread y be en experimen- 
tally realized |H [H Oil . In Ref. |2fl, a vortex line 
with one unit of circulation was imprinted on a pair of 
condensates occupying different internal levels of 87 Rb 
using a Raman coupling which was set locally resonant 
in a small region by rapidly rotating laser beams. In 
Refs. |27l |23, |29j , singular singly-quantized and doubly- 
quantized, as well as non-singular coreless vortex lines 
were created by adiabatically inverting the magnetic bias 
field along the trap axis. Also dark solitons have been im- 
printed on atomic BECs, e.g., by imaging the atom cloud 
through an absorption plate [20, EI ■ Such n phase kink 
planes have been observed to decay into a hybrid of vor- 
tex rings and lines through the dynamical 'snake' insta- 
bility nmm. a controlled method of creating vortex 
rings and particlelike solitons was proposed in Ref. 0, 
where it was shown that a topological phase singular- 
ity forming a closed circular loop may be imprinted on 
the matter field while changing the internal state of the 
atoms. Ring defects in spinor BECs may also form from 
simpler core structures as a result of dissipation 33] . The 
defect generation using a 'light roadblock' was first real- 
ized in Ref. 0. More complicated defect configurations 
were recently produced by superposing multiple road- 
blocks [H|]. It has also been proposed that the strong 
light-matter coupling of the ultra-slow and stopped light 



[3~i. I35L l3l| could be used to efficiently transfer vortex 
states between light and atomic BECs and to store light 
modes with orbital angular momentum inside the BECs 
|2*H , as well as to use a vortex lattice to create a photonic 
band gap |37|. 

In Section [n] we show how em fields containing phase 
singularities can be used to imprint vortex rings on BECs. 
We discuss several examples of how the required fields can 
be generated by superpositions of plane and/or Gaussian 
laser fields. In Section UTTI we apply the method to gener- 
ating two vortex rings and studying their ensuing recon- 
nection dynamics. In Section llVI we then turn using the 
light roadblock to generate of systems of multiple vortex 
point particles in a 2D geometry, focusing on collisions 
of vortex pairs, which dissipate energy via phonons and 
sometimes leads to vortex annihilations. 



II. CONTROLLED PREPARATION OF 
VORTEX RINGS 

A. Coupled two-component condensate 

A vortex ring can be engineered by using an em field to 
imprint topological phase singularities on the matter field 
while changing the internal state of the atoms |3| . Here 
we consider a two-component BEC where the two inter- 
nal states are coupled by means of the em fields with the 
Rabi frequency, fi(r). The dynamics of the BECs with 
the Rabi coupling between the levels \i) and \j) follows 
from the coupled Gross-Pitaevskii equation (CPE) 

ih^ = (H Q + Si +J2 Kik\ipk\ 2 )ipi + fin*ipj ■ (1) 
k 

We assume the two-component condensate, with the to- 
tal number of N atoms and the atomic mass m, to be 
confined in a perfectly overlapping, isotropic trap with 
the trap frequency lo: 

h 2 I 

H = -—V 2 + -mwV . (2) 
Zm I 

The parameter Sj(r) incorporates the detuning of the 
em fields from the resonance of the internal transition 
as well as em field-induced level shifts which may be 
non-uniform in the case of a spatially varying field in- 
tensity. We have also defined the interaction coefficients 
Kij = 4irh 2 aijN/m, with the intraspecies and the inter 
species scattering lengths denoted by an and a,j (i =/= j), 
respectively. Such a system has been experimentally re- 
alized using, e.g., the |2) = \S\/ 2 ,F = 2, Mp = +1) and 
|1) = \S\/2,F = 1,Mf = —1) hyperfine spin states of 
87 Rb. The two hyperfine levels were coupled by means of 
a two-photon transition (one microwave and one rf pho- 
ton). For the 87 Rb components the interaction strengths 
are nearly equal, with an : a 21 : a 22 1-03 : 1 : 0.97 and 
a 21 ~ 5.50 nm 38]. Since the scattering lengths satisfy 
a 2i ~ a n a 22, the two species experience dynamical phase 
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separation and can strongly repel each other [38j . In the 
absence of the em coupling between the different compo- 
nents the interatomic interactions of the two 87 Rb com- 
ponents do not mix the atom population and the atom 
numbers of the both species are separately conserved. 

The phase profile of the driving em field can be im- 
printed on the matter wave by means of transferring the 
atomic population between the two different internal lev- 
els. We assume that all the BEC atoms initially occupy 
one of the hyperfme states, let's say |2). Some popu- 
lation is then transferred from \2) to |1) by means of a 
Rabi pulse f2(r, t). After the pulse, the relative phase be- 
tween the BECs in the two levels is proportional to the 
phase profile of the Rabi field, as indicated by Eq. @J. 
We may construct a phase profile for the em field, where 
the node points of the field amplitude may correspond to 
the topological singularities of the phase of the em field. 
The em coupling provides then a method for imprinting 
these topological singularities on the condensate. In or- 
der to imprint a vortex ring on the z — plane, centered 
at the z axis, we require the em field to exhibit a closed 
circular phase singularity in such a way that in the close 
neighborhood of the circular node the field is of the form 



fl(r, t) ~ K(t) [(p - pa) + ijz] = n(t) r)e 



it) 



(3) 



where p = (x 2 + y 2 ) 1 ^ 2 and po denotes the radius of the 
ring. Here we have also defined 7/ = [(p — po) 2 + 7 2 z 2 ] 1/ ' 2 
and 9 = arctan [yz/(p — po)]- The field amplitude © 
vanishes at the ring (p = po, z = 0) with the desired 2tt 
phase winding along any closed loop encircling the ring, 
representing one unit of quantized circulation around the 
ring. The parameter 7 describes the anisotropy of the 
vortex core. The value 7=1 corresponds to an isotropic 
core. 



B. Imprinting a vortex ring 

A technique to imprint a vortex ring on the BEC was 
proposed in Ref. by means of constructing an em field 
amplitude which exhibits the desired form ©. This was 
obtained using 



Q(r) = fi [a-f(x,y)+i/3sm(kz)] 



(4) 



Here the Rabi amplitude represents a coherent superpo- 
sition of a standing wave along the z axis, a constant 
field a, and the field f(x,y), which could be a Gaussian 
field focused weakly on the xy plane or, alternatively, a 
superposition of two standing waves along the x and the 
y directions: 



f(x,y) = exp - 



f(x,y) 



£1 

e 

2x 

2\ C0S T 



cos ■ 



2y 



(5) 
(6) 



We assume p/£ ~ \J\ — a <C 1. Then by expanding 
Eqs. 0)-© first order in p/£ and k\z\ <C 1, we obtain 
Eq. with po = — o- and 7 = /3£ 2 fc/2p . Hence, a 
constant field with three orthogonal standing waves, or, 
alternatively, with a standing wave and a parallel Gaus- 
sian beam, are sufficient to imprint a vortex ring on an 
atomic BEC. Moreover, by choosing a = 1 — (/3£/c/2) 2 , 
yields an isotropic vortex core with 7=1. 

The field configuration Q is most suitable for mi- 
crowave (or longer wavelength) em fields for which the 
wavelength A = 2n/k is longer than the typical radius 
R of the BEC. For optical fields A < R the em cou- 
pling would create multiple copies of the ring, displaced 
from each other by A. However, even for optical fields it 
should be possible to shape the wave fronts of the cou- 
pling lasers in order to avoid rapid phase variation at 
the length scale A. This could be done, e.g., by using 
microlens arrays |39l | , laser beams co-propagating in the 
direction of the standing wave z, or beams with the wave 
vectors nearly perpendicular to z. Here we show that 
another alternative is to consider two-photon transitions 
via some intermediate level |3) and use slightly different 
wave numbers for the lasers. Consider the first transi- 
tion 1 2) — > 1 3) to be induced by four co-propagating laser 
beams, one of which is weakly focused on the xy plane, 



r> 2 (r) =^2 



a — e p e 



A f p i(k+k')z _ A{k-k'); 



(7) 



Here we assume that k' <C k and the coupling frequency 
□2 = d23 • E ^ /h is in this case determined in terms of 
the atomic dipole matrix element d23 of the transition 
|2) — > |3) and the positive frequency component of the 
driving electric field Ej. The second transition |3) — > |1) 
is driven by a non- focused co- propagating field f2i(r) = 
f2i exp (ikz), which cancels the rapid phase variation of 
along the z axis. In the limit of large detuning A 
of the laser $12 (r) from the resonance of the |2) — > |3) 
transition, the effective Rabi frequency for the |2) — > |1) 
transition then reads 



n 2l (r) 



2n 2 (r)nj(r) 
A 

2n 2 ni 



+ i/3 sinfc'z 



(8) 



This has the desired form In addition, even though 
we may have kR » 1, it should be possible to choose 
k' k, so that k'R < 1. This allows us to create pre- 
cisely one vortex ring in the BEC. 

In a more challenging scheme, one may construct an 
optical lattice with a long periodicity [l!| by means of 
superposing two separate two-photon transitions 



sin[(fci - fc 2 ) 



= — (e lklZ e 
2i v 



(9) 



such that (ki — k 2 )R < 1. Here each exponent factor 
represents a running wave transition by one photon. 



In the following, we demonstrate how the vortex ring 
could also be imprinted using only two Gaussian laser 
beams providing an experimentally simple scheme for a 
controlled creation of vortex rings. 



C. Imprinting a vortex ring using Gaussian beams 

Here we show that a vortex ring defect can also be im- 
printed on an atomic BEC in a controlled way by means 
of two parallel Gaussian laser beams. The advantage of 
this scheme, as compared, e.g., to Eq. Q is the obvi- 
ous simplicity of the field configuration and the common 
experimental availability of the beams. The complete 
expression of the Gaussian laser beam reads: 



G{r,w Q ) 



_ WQ 
w(z) 

x exp 



exp 



ikz 
2 i 



z .. p 

i arctan h ik 

z 2R(z) 



w(z) 2 



(10) 



where w(z) = tUo"\A + z2 / z o an< ^ R( z ) = z (l + z o/ z2 )- 
Here wo represents the minimum beam waist and zq = 
ttwq/X is the Rayleigh range, the distance over which the 
focusing is increased to \/2wo due to diffraction. In a 
weakly focused beam, for p <C Wo, \z\ <C zq, Q{t,Wq) ~ 
exp {ikz — p 2 /wq). 

We take the two Gaussian beams as a phase-coherent 
superposition: 



Q r (r) = Q [<?(r, w ) + cQ(r, w' Q )] , 



(11) 











c = — exp 


-pl 







for some complex c. It is straightforward to see that by 
choosing 



(12) 



the field configuration (|ll(l exhibits a circular node at 
(p = po, z = 0) with the 2ir phase winding around the 
ring. The complete destructive interference of the two 
Gaussian beams does not occur outside the plane of the 
minimum beam focus size due to diffraction. By modify- 
ing the relative amplitude and the focusing of the beams, 
the size of the ring can be changed in a controlled way. 
The anisotropy of the core is determined by 



WqWq 2 - K +Wq)p 2 o 
w^w' 2 p k 



(13) 



In order to avoid the rapid phase variation of the 
Gaussian laser along the z axis to be imprinted on 
the BEC, it is again advantageous to consider a two- 
photon transition via some intermediate level |3) using 
co-propagating lasers, so that, e.g., the field O r (r) in 
Eq. (|llfl drives the transition |2) — ► |3) and a non- focused 
co-propagating field O31 exp (ikz) drives the transition 
1 3) — > |1). In the limit of large detuning A, the effec- 
tive Rabi frequency for the |2) — > |1) transition then 




FIG. 1: The amplitude of the superposition of two Gaussian 
laser beams at z — (on the left). The cylindrically symmet- 
ric amplitude exhibits a node point at po — 1.51. The constant 
surface plot of a vortex ring in an atomic BEC (on the right) 
with the phase profile imprinted from the two Gaussian laser 
beams. The Rabi frequency tx e~ lk2Z Q r . Here the beam pa- 
rameters are wo = 3.0/, w' = 4.5Z, kl = 7, k-^l = 6.86, and 
the radius of the ring po = 1-51. As a result of the imprinted 
phase pattern, the vortex ring slowly moves along the beam 
axis. 



reads 2£l 31 Q r (r) exp (— ikz)/ A. The corresponding spa- 
tially non-uniform laser-induced level shifts are <5i(r) ~ 
2|fi 31 (r)| 2 /A + (5 12 and S 2 (r) ~ 2|ft r (r)| 2 /A, where 5 12 is 
the effective two-photon detuning from the resonance of 
the internal transition |2) — > |1). 

In Fig. n we show an example of a vortex ring in a 
two-component BEC, obtained by imprinting the phase 
profile of the two Gaussian laser beams from Eq. (|11|) . We 
use the parameters of 87 Rb with the nonlinearity K21 = 
b OOhuil 3 and the radius of the ring po = 1.51, where I = 
yjh/muj. The ring occupies the level |1) with the atom 
numbers N2/N1 ~ 6.2. For the chosen parameters, the 
vortex core is noticeably anisotropic with 7 ~ 0.06. A 
more isotropic core may be obtained, e.g., by decreasing 
the radius of the ring by means of changing the value c 
according to Eq. (|12|l . For po = 1 .0Z we obtain 7 ~ 0.12 
and for po = 0.5/, 7 ~ 0.27. As a result of the nonlinear 
evolution and dissipation, an anisotropic vortex core in 
an atomic BEC is generally expected to relax towards an 
isotropic core shape. 

Experimentally, one could prepare the phase-coherent 
field configuration (| 1 II) using a single laser source which 
is split in two and later recombined after appropriately 
modifying the relative amplitude and the beam focus- 
ing. One possible experimental limitation might then be 
a slight misalignment of the two recombined Gaussian 
beams. We also tried field configurations with the prop- 
agating axes of the two beams displaced. It is quite easy 
to show that the vortex ring singularity Eq. I jlljl is robust 
against perturbations where the two Gaussian beams are 
slightly modified. 



D. Multi-beam superpositions 

Two beam configurations, such as Eq. (|llf) . or more 
complicated multi-beam superpositions can also be pre- 
pared using diffractive optical components. In particular, 
computer-generated holograms and spatial light modu- 
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lators may be used to prepare the desired optical field 
superpositions which are required to imprint the vortices 
on the atomic BECs. By using a third Gaussian beam 
in Eq. (|llf) we may better optimize the field parameters 
according to a steady-state vortex ring solution, or, al- 
ternatively, to engineer some desired excitations in the 
ring. With the superposition: 

J2 r (r) = Q {G{r,w ) + Cl g{r, Wl ) + c 2 g(r,w 2 )] , (14) 

we may, e.g., choose the parameters c\ and c 2 according 
to the given vortex ring radius and the core isotropy, 
while wo, wi, and w 2 may be optimized to determine the 
spatial profile of the BEC wave function. 



E. Knotted vortex lines 

More complicated field superpositions than the one 
shown in Eq. fTljl have been successfully created using 
spatial light modulators. In Refs. [if! Elj a field con- 
figuration of four optical Laguerre-Gaussian beams, with 
each exhibiting zero, one, or two units of orbital angular 
momentum, were generated using a liquid crystal array 
acting as a phase mask and controlling both the field 
amplitude and the phase. Such a field superposition pos- 
sesses a phase singularity forming a closed knotted loop 
|42t l43l | (a torus knot) in a paraxial field and, in princi- 
ple, could also be imprinted inside a BEC, provided that 
the hologram has a sufficient spatial resolution in order 
to allow the beams to be focused inside the atomic cloud. 



III. ENGINEERING VORTEX 
RECONNECTIONS 

A. Preparing multiple vortex rings 

We may use the techniques described in the previous 
section also to imprint multiple vortex rings on the BEC. 
Generally this can be done by using multi-photon tran- 
sitions, where the effective Rabi pulse is of the form: 

n 

fi(r, t) ~ K(t) ni(P (j) - P ( o } ) + - 4 j) )] Pi > (15) 

where we write p& = [(x^ ~ x[ j) f + (j/W - y^) 2 ] 1 / 2 , 
representing n vortex rings on the plane z^ = 0, cen- 
tered at (x^ = xtf\y® = y { o\z^ = z { j) ), with the 
radius given by p^\ All the vortex rings in Eq. I|15|) 
do not need to have the same orientation, but the axis 
r"' may involve different spatial rotations with respect 
to some fixed axis r. The exponent pj in Eq. I|15[) denotes 
the topological charge of the vortex ring j. 

We numerically studied the preparation of two 
nonoverlapping vortex rings by integrating the coupled 
GPE ||TJ in the presence of the em coupling. In Fig. [21 



we show an example of two vortex rings with the radii 
pa = 0.91 prepared on the z = plane. Both rings are 
displaced from the trap center in opposite directions by 
1.5Z. Here the vortex rings are created by means of the 
em field f2(r — ri)0(r — r 2 ) (for ri = — r 2 = 1.5ly) with 
fi(r) determined by Eqs. (0J and (JSJ). The wavelength 
A = 261, (3 ~ 8.3 x 1CT 3 , and the width £ = 30Z. We use 
the parameters of 87 Rb with the atoms initially occupy- 
ing the level |2). The nonlinearity is K21 = 430Jtu1 3 . A 
short pulse tu> = 0.1 is applied to transfer population to 
the level |1) in order to prepare the vortex rings. This 
could be, e.g., a two-photon transition via some inter- 
mediate level 1 3), where the field f2(r — ri) drives the 
transition \2) — > |3) and the field f2* (r — r 2 ) drives the 
transition |3) — * |1). Alternatively, we could have used 
the Gaussian laser beams ()llf) . e.g., driving a two-photon 
A three-level transition, where the beams inducing the 
second transition are co-propagating with the first pair 
of beams, but are multiplied by the phase factor (—1) and 
displaced in space by the amount of 3Z. By modifying the 
strength and the duration of the Rabi pulse we can con- 
trol the excitation of the vortex ring pair. This allows us 
to study the dynamics of the vortex ring interactions for 
different initial conditions. 



B. Vortex ring reconnections 

Vortex reconnections represent basic vortex-vortex in- 
teractions. These are of fundamental interest in their own 
right, as well as providing a useful model to investigate 
superfluid turbulence which can be considered as being 
due to the dynamics of a dense tangle of quantized vor- 
tex filaments . The evolution of such a tangle crucially 
depends on the vortex reconnections. The vortex recon- 
nections also have an analogy with the models describing 
the evolution of cosmic strings in the early Universe cos- 
mology At very low temperatures, with negligible 
normal fluid fraction, the superfluid turbulence and the 
vorticity is expected to decay as result of sound emission 
|45| . Sound emission may occur due to vortex motion or 
reconnection. 

The proposed techniques of vortex imprinting provide 
a useful tool in investigating vortex reconnections in su- 
pcrfluids. Atomic BECs can be cooled down to very 
low temperatures with negligible normal fluid compo- 
nent. By modifying the pulse parameters of the imprint- 
ing process we may also control the level of excitation in 
the system. Moreover, vortices in two-component BECs 
offer several additional advantages. First, it is possible 
to monitor the system in real time: The filling of the 
vortex cores by the other BEC significantly increases the 
core size, as compared to vortices in a single-component 
BEC, and makes them observable in situ, even without a 
ballistic expansion. Second, the nonlinearity of the vor- 
tex excitations can be easily controlled by varying the 
relative density in the two components (via the strength 
and the duration of the coupling laser fields). 
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The vortex reconnections have previously been mi- 




FIG. 2: The snapshot images of the isosurface plots during 
the time evolution of two merging vortex rings in a trapped 
two-component BEC. We show the BEC component contain- 
ing the two vortex rings, while the other component exhibits 
a uniform phase profile. The initial state was prepared by 
means of an em coupling field between the two BEC compo- 
nents, which imprints the desired phase singularities on the 
matter field while changing the internal level of the atoms. In 
order to relax the energy of the initial state, we integrated the 
wave functions in imaginary time for the duration of tuj — 0.06 
before the real time evolution. We show the isosurface plots 
at times tu = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 1.0, 1.4. 



EH l4p in a homogeneous single-component BEC. In 
Refs. |48ll49l | two vortex rings with an initial relative ve- 
locity were scattered from each other, varying the angle 
of intersection and the impact parameter. As a result, a 
substantial loss of vortex line length was observed, which 
was attributed to sound emission. 

As an application of the multi-ring imprinting tech- 
niques, we showed here how to create two initially static 
vortex rings in a trapped inhomogeneous two-component 
BEC, which will later interact and collide. The two vor- 
tex rings on the same z = plane were imprinted using 
cm fields, as explained in the previous section. Due to the 
filling of the vortex core by the second BEC component, 
an isolated vortex ring would approximately remain in 
the same location in the trap on the time scale of the 
trap period. If the initial separation of the rings is suffi- 
ciently small, the vortex rings are attracted to each other 
and merge locally, after the em coupling fields have been 
turned off. 

Qualitatively, the long-range interaction energy be- 
tween vortex filaments s and s' with the vorticity k and 
k! may be obtained, analogously to the magnetostatics, 
from [5(| 

tikk' f ds-ds' 
E -^TJ \r(s)-r( S i)\- (16) 

Two antiparallel vortex filaments are attracted to each 
other. Once the vortex cores are close enough, they 
merge at the point of the closest approach. Since the 
facing segments of the two similarly circulating rings on 
the same plane represent vorticity with the opposite cir- 
culations, these segments are attracted to each other and 
eventually they annihilate. Consequently, the two rings 
are joined together. By integrating Eq. (jlbfl with an ap- 
propriate cut-off, we may also obtain the energy of an 
isolated vortex ring as E cx npo[lri(epo) — 2], where n de- 
notes the density and the parameter e depends, e.g., on 
the core thickness |43, E3 and where po denotes the ring 
radius. Close to equilibrium the reduction of vortex core 
length can indicate the emitted sound energy |52j . 

In the numerical studies of the two-component BEC 
dynamics, we use initial states corresponding to different 
strength of excitation. These could be created by chang- 
ing the strength and the duration of the em pulse in the 
imprinting process. However, for computational simplic- 
ity, here we always use the same pulse parameters, but 
relax the initial state by means of varying length of imag- 
inary time evolution before the actual dynamics. 

In Figs. [21 and [3] we show the snapshot images of the 
real time evolution of the two vortex rings after they were 
prepared by means of the Rabi coupling, as explained 
in the previous section. The images were obtained by 
numerically solving the coupled CPE for the trapped 
two-component BEC in the absence of the coupling fields 
ft. The integration was performed using the split-step 
method [5i| on a spatial grid of 128 3 . The 3D mean- 
field dynamics corresponds to the parameters of 87 Rb 
experiments, as described in the previous section. The 
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number of atoms in the two components N% — 0.867V 
and JVi ~ 0.147V, where the total atom number N — 
430V(47rai 2 ). 

In Fig.|21the initial state, obtained from the numerical 
integration of the Rabi coupling, was relaxed by means 
of imaginary time evolution of duration tu = 0.06 before 
the displayed real time evolution. In Fig. [3] the initial 
imaginary time evolution was t<u = 0.02 and the merging 
process of the two rings is more violent. During the imag- 




FIG. 3: The same system as in Fig.|3 but the initial state of 
the real time evolution was obtained by means of propagating 
the wave functions in imaginary time for tu) — 0.02. The 
resulting reconnection process is more violent than the one 
displayed in Fig. [5] We show the isosurface plots at times 
tu = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0. 



inary time evolution we separately normalized the both 
wave functions tpi and ^2 > corresponding to the separate 
conservation of the atom number in each BEC compo- 
nent in the absence of the coupling field. We also show 
the total condensate energy density et , averaged over the 
spherical angles, 

= I dyd{cos6) £ iR(H + ^|V>i| 2 )V>i , 

(17) 

as a function of the BEC radius and evolution time in 
Fig. 0] for the cases shown in Figs. |2] and |3J as well as for 
cases without the initial imaginary time evolution and for 
the imaginary time propagation of duration tu> = 0.1. In 
Fig. [S] we show the contribution e\ to the energy density 
which only depends on the atom density in the compo- 
nent |1) containing the vortex rings: 

£l (r) = i- j d<Mcos0)Vh*(#o + ^|^i| 2 )V>i- (18) 

In Figs. [5] and [3] the vortex rings are attracted to each 
other as their motion is also bent downwards. It should 
be emphasized that in the two-component BEC the vor- 
tex rings are initially static, there is no thermal atom 
component, and, after the preparation of the initial state, 
the dynamics is purely Hamiltonian with conserved total 
energy and no added dissipation term. In order to vortex 
rings to merge, they need to lose energy. In Figs. 14151 we 
can clearly identify the released excess energy as sound, 
radiated outwards and transformed into surface excita- 
tions. It is the turbulent dynamics triggered by the re- 
connection process which locally produce an effective dis- 
sipation mechanism for the vortex configuration to relax. 
We can enhance the sound emission and the reconnec- 
tion dynamics by considering a more rapid imprinting 
process. It is interesting to compare the merging vortex 
rings to the formation of a vortex lattice in a rotating 
BEC. In the latter case turbulent dynamics is triggered 
by the rotating potential and the vortex lattice can lo- 
cally relax in the classical mean-field dynamics, even at 
T = and without any explicit damping term 0. 

The reconnection process in Fig. 0] begins with a sig- 
nificant concentration of energy towards the trap cen- 
ter, as energy is propagated towards r = 0. The vor- 
tex motion and the vortex reconnection process excite 
Kelvin modes, resulting in a gradually decreasing vortex 
line length. The Kelvin mode excitations of the vortex 
cores are clearly observable in the first snapshot images 
of Figs. 12131 The generated vortex cusp at the reconnec- 
tion region detaches from the merging vortex rings and 
is emitted as sound radiation. In the process the vortex 
energy is converted into a sound pulse which propagates 
away from the vortex. We observe large emission of vor- 
tex rings and/or rarefaction waves [54( at the reconnec- 
tion region, e.g., in Fig. [3at time tu = 0.6 and in Fig. [21 
around 0.5 < tu < 0.6. The emission events can also be 
seen in Fig. [SJas high energy regions close to the origin. 
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FIG. 4: The BEC energy density et(r,t) averaged over the 
spherical angles during the reconnection process. The initial 
imaginary time evolution cot — (top left), uot = 0.02 (top 
right), Lot — 0.06 (bottom left), and Lot — 0.1 (bottom right). 
The white color indicates a high value of the energy density. 
The sound waves first propagate towards the trap center when 
the vortex cusp is generated as the two rings merge. The 
energy is then radiated outwards and is transformed to surface 
excitations. The straight vertical stripes are an artifact of the 
numerical averaging procedure over the spherical angles. 
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FIG. 5: The energy density ei (r, i) averaged over the spherical 
angles, depending only on the density of the BEC component 
containing the vortex ring. We show the same cases as in 
Fig. 2] In the second and the third figure one may clearly 
recognize the energy of the emitted sound pulses in the fifth 
and the sixth images of Figs. |2] and |3] 



The emitted sound energy results in strong surface ex- 
citations of the BEC (displayed in the last two images 
of Figs. |21 and |2l with the fastest sound pulses (shortest 
wavelength phonons) reaching the surface first. 



IV. USING ULTRA-SLOW LIGHT IN THE 
STUDIES OF VORTEX COLLISIONS 

A. Dissipation in vortex pair collisions 

One of the fascinating aspects of vortices in current 
BECs is the fact that the relatively small nonlinearity 
leads to many finite compressible effects. In geometries 
where the system is essentially two dimensional, due to 
tight confinement in one direction (we will label y here), 
vortices act as point particles. In an incompressible fluid, 
they are indestructible particles which interact via con- 
servative potentials, moving always according to the ve- 
locity fields created by neighboring vortices. For exam- 
ple, two vortices of like charge will orbit each other at 
some constant radius, while oppositely charged vortices 
will form 'pairs' which travel along straight lines main- 
taining the same distance between them d |15| . In ei- 
ther case, the velocities are inversely proportional to the 
distance to the neighboring vortices. However, when the 
distances between vortex pairs become comparable to the 
healing length (which is the characteristic size of vortex 
cores), this picture breaks down. For example, in the case 
of co-rotating vortices of like charge, the pair will emit 
sound radiation and slowly drift to larger distances, low- 
ering their energy |55| . The radiation occurs due to accel- 
eration of the vortices as they undergo circular motion. 
Finite compressibility effects have also been important in 
explaining discrepancies between experimental and the- 
oretical descriptions in excitations of vortex lattices [56( . 

The dynamics of multiple vortex configurations, partic- 
ularly ones containing vortices of both circulations, can 
be particularly rich. For exam ple , scattering between 
vortex pairs have been studied (53 as well as the exis- 
tence and behavior of soliton-like structures in 2D geome- 
tries, and how they can dynamically form from closely 
spaced vortex pairs |5^|. However, there has been little 
work on cases where the dissipation due to compressibil- 
ity is significant and vortex annihilations occur. The dis- 
sipation is analogous to the length shortening and sound 
emission which occurs with colliding rings in 3D geome- 
tries studied above and in Rcfs. 

mm 

Here we show that the recently developed technique of 
the light roadblock, based on ultra-slow light propaga- 
tion in atomic BECs, may provide a method of prepar- 
ing point vortices and systematically investigating their 
collisions in strongly anisotropic 'pancake-shaped' atom 
clouds. We review the ultra-slow light propagation in 
BECs and the creation technique by which systems con- 
taining multiple vortices of each charge can be engi- 
neered, allowing study of these interesting phenomena 
which have remained unexplored experimentally to now. 
In particular, we show how it is possible to induce vortex 
pair scattering, in a regime where the dissipation due to 
finite compressibility is quite apparent. We then study 
such scattering in more detail to obtain new quantitative 
results on the role of dissipation in such collisions. 
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B. Ultra-slow light propagation and the light 
roadblock 

Slow light is based on electromagnetically-induced 
transparency (EIT) [Toj . a quantum interference effect 
that permits the propagation of light through an other- 
wise opaque medium. A coupling laser, resonant with 
some initially unoccupied hyperfine state |2) and an op- 
tically excited state |3), is used to create the interfer- 
ence resulting in transparency (vanishing resonant ab- 
sorption), a very large dispersion, and ultraslow group 
velocity of a probe pulse, resonant with an initially oc- 
cupied state |1) and |3). Probe light pulses have been 
slowed down and spatially compressed (compared to their 
free space propagation values) by up to eight orders of 
magnitude 0,H,l3- As tne probe propagates through 
the BEC, it puts the atoms into spatially dependent su- 
perpositions of |1) and 1 2) with the relative phase and 
amplitude of the wavefunctions (the atomic coherence) 
reflecting the amplitude and phase pattern of the input 
probe. 

Combining this technique in BECs with the recently 
developed experimental technique of the 'light roadblock' 
, whereby a fast spatial variation of the coupling field 
is introduced, provides an efficient method of inducing 
phase singularities on matter waves. Sharp density de- 
fects are prepared by abruptly perturbing the light prop- 
agation. In anisotropic 2D trap configurations the cre- 
ated defects evolve in a predictable way to solitonlike 
phase kinks and subsequently then to systems of multi- 
ple point vortices (of both circulations) whose collisions 
may be investigated in detail. Here we focus on the case 
of a strongly anisoptropic trap, leading to essentially 2D 
dynamics. In Ref. 6], the geometry was truly 3D and 
in Ref. [33 the 3D aspects of the resulting topological 
structures in a similar system were investigated. 

We specifically consider the BEC occupying two hyper- 
fine levels |1) and |2). These are connected, respectively, 
to a common excited state |3) by orthogonally propagat- 
ing, resonant laser fields: the +z propagating probe f2i3 
and +y propagating coupling fl23, each of wavelength A. 
The excited state |3) decays at T, forming a A three- level 
structure. In the slowly varying envelope (SVE) approxi- 
mation for the light fields Cta (with the rapid phase rota- 
tion at the optical frequencies and optical wavenumbers 
factored out) the Maxwell's equations read |Hflj : 

+ n 23 r^2e- lk( - z - y) ], (19) 

+ n 13 r 2 ^ie ik{z - y) }. (20) 

Here N is the initial total number of BEC atoms, k = 
27r/A, /,3 are dimensionless oscillator strengths, and a = 



3A 2 /27r is the resonant cross-section. In Eqs. I|19l20ll we 
have adiabatically eliminated tp 3 by assuming that the 
dynamics of the internal atomic degrees of freedom and 
the light are much faster than the external dynamics |59l 
l60j | . The BEC wavefunctions tpi , 1^2 evolve according to 
generalized GPEs: 

ihtpi =(h + ^ K ik\ipk\ 2 )i>i 

k 

- [|a| 2 ^ + njfyVie-**^-^] , (21) 

where i = 1, 2, j ^ i and ej = z, ti — V- The last term 
in Eq. (|21|l results in both coherent exchange between 
|1), |2) as well as absorption into |3). The phase factors 
assure that atoms coherently coupled into |2) experience 
a two-photon momentum recoil relative to the nearly sta- 
tionary |1) component. In our model, atoms which pop- 
ulate 1 3) and then spontaneously emit are assumed to be 
lost from the BECs. 

If one were to input a weak (|Oi3 1 <C |^23|) free- 
space pulse of some length L z into a BEC, simultane- 
ous solution of Eas. (|19l21fl reveals that the pulse is com- 
pressed to a length (V g /c)L z , where the group velocity 
is V g — |f223| 2 /r/i3t7iV|'i/>i| 2 , as it propagates through 
the BEC. Any phase and amplitude features along the 
longitudinal direction are similarly compressed while the 
transverse features are unaffected. While the usual de- 
scription of slow light propagation depends on the weak 
probe assumption ( | ^2 13 1 <C 1^23 1), the basic transfer of 
amplitude and phase patterns works even when this as- 
sumption is not strictly satisfied. This regime was stud- 
ied in detail in Ref. |59|. 

Spatial modulation of the input coupling field along z 
is then used to vary the speed and length of slow light 
pulses as they propagates through the BEC. In Ref. 6] 
this variation was accomplished by using a razor blade 
to block the coupling field from the BEC in the z > 
region so near z — the coupling intensity quickly varied 
from some large value to zero. In this region, the probe 
beam is slowed to zero group velocity and the length be- 
comes arbitrarily small. At the same time the probe to 
coupling intensity ratio becomes large in this region, and 
so the fraction of atoms coupled into |2) becomes large. 
The result is a very narrow, large amplitude, density de- 
fect in the original condensate internal state |1), with the 
atoms in this region coupled into the other stable inter- 
nal state 1 2). Because of the large photon recoil in the 
orthogonal geometry, atoms transferred to |2) leave the 
defect region in < 1 ms, leaving an empty hole in a sin- 
gle component BEC. This method is particularly useful 
for inducing density features in BECs smaller than the 
experimental optical resolution of the system, as the spa- 
tial compression at the roadblock occurs due to the slow 
light propagation effects and not focusing. This is partic- 
ularly valuable in creating topological defects, since their 
characteristic length scale, the healing length, is typically 
< 1 //in. Numerical simulation of the 23 Na experiment in 
Ref. [|| indicates that density defects ~ 2 \xm long (and 
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nearly 100% depth) were being created even though the 
spatial variation of the coupling field (which is limited by 
optical resolution) was about ~ 10 fim. This interpreta- 
tion was further supported by the fact that the density 
defects of the features were unobservable in in-trap ab- 
sorption images but were clearly visible after 1 ms of free 
expansion. 



C. Inducing vortex pair collisions with a 'light 
roadblock' 

In the original experiment the imprinted density 
defects were seen to subsequently break up into a series 
of a solitons, due to the superfluid analog of a shock 
waves. These solitons, in turn, broke u p in to pairs of 
vortices via the snake instability [nj, ll2L Il3j] . We show 
an example of how this comes about in Fig. 

In this example we consider a BEC with repulsive 
interactions in a trap with much stronger confinement 
along the y axis, with the trap frequencies satisfying 
u) y 3> uo x , uj z , so that the 2D limit of the GPE is valid. We 
assume that the defect has been created using the light 
roadblock technique, as explained earlier. Moreover, af- 
ter the light fields have been turned off, we assume that 
the BEC component |2), which is initially filling the den- 
sity defect in |1), is instantaneously coupled out of the 
trap due to the photon-induced recoil, as in the origi- 
nal experiments 6]. Then the BEC dynamics approxi- 
mately follows from the single-component 2D GPE for 
ip(x, z), when we write ^i(r) ~ (j){y)ip(x, z), in terms of 
the ground state harmonic oscillator wave function <f>(y) 
along the y axis: 

t^ = [~^+^^+^) + ^u>\^, (22) 
ot 2m 2 

where k 2 e> — Kii / I y y2/K , with l y = yjh/muy. The 
additional density-dependent contribution to the scat- 
tering length in 2D is negligible when \/2~%l v /a ^> 
In (8n 3 / 2 l y n 2D a) H2], where n 2D = N]^] 2 denotes the 
2D density. In the numerics, we choose uj x /oj z ~ 3.81 
and the nonlinearity n 2 £> = 2360 hui z l 2 (l z = \J h/muj z ). 
The numerical results depend on the dimensionless non- 
linearity R 2 d = n 2 D/huJl = 2\/2~jTNa/l y , so our results 
are unchanged in any rescaling of the parameters which 
maintains constant Na/l y and l z /l x , provided that the y 
confinement remains large enough for the 2D GPE i|22|) 
to be valid. 

We then remove atoms from the relaxed ground state 
to create a narrow (~ 0.7/ z half width) density hole, with 
100% of the density removed along the central z = axis. 
The GPE ill' 21) was propagated with a Crank-Nicholson 
grid with an equal spaced grid (300 points in x and 600 
in z), with time steps 1.65 x 10~ 7 /w z . Both were varied 
to assure the results were not effected. 

To make the connection with experimental parame- 
ters, we quote our results in terms of 23 Na parameters, 



-40 -20 20 40 -40 -20 20 40 



1 1 

-0.0 ^ 

i i 1 




1 1 1 1 1 

i i i i 


l l | 

- ^Ml 

5.0 ms 

i i i 


i i 

m ; 


i i [ i i 




i i i 

_i5.orai^*"' 

. i i 


► ; 


II 11 

_79.o" mP^* 




' i 

i 


* : 


i 


' ' 1 - ' ' 


i i | i i 









I I I I I I I I I I I 

-40 -20 20 40 -40 -20 20 40 
z(nm) z (nm) 

Dens. (10 4 cm' 2 ) q ^^^H 8.6 



FIG. 6: The evolution of a 3 fim (1/e half-width), full depth 
defect in a sodium BEC at the times indicated. Qualitative 
features of the evolution are explained in the text. The density 
scale refers to the normalized 2D density \ip(x, z)\ 2 . 

which has a scattering length a ~ 2.75 nm, and choose 
a trap oj z = 2ix x 21 Hz, giving l z ~ 4.6 /xm and de- 
fect size ~ 3 fim. The ground state BEC in this trap 
has a central peak density of n^n = 1-5 x 10~ 14 cm" 3 
and the superfluid healing length at the cloud center is 
C ~ h/{2mn 2D \^\ 2 ) 1 / 2 ~ 0.4 fim. 

In the limit of small-amplitude and large wavelength, 
the imprinted defect would split into two counter- 
propagating sound waves [6^J . However, due to the non- 
linearity of the excitation (and the small l engt h scale) @ 
it here sheds off a series of a gray solitons 63 , as seen in 
the snapshot at 5 ms. In addition, the slowest, and deep- 
est soliton has begun to curl due to the snake instability, 
and vortices have begun to form. Because angular mo- 
mentum is not being imparted to the system, vortices are 
formed in pairs of opposite charge. Investigation of the 
phase shows all the vortices formed are singly quantized. 

At 15 ms each quadrant contains a line of three vor- 
tices (alternating in sign along each line) . Note that they 
are closely spaced here so their density profiles overlap, 
and the identification of the vortices is done by finding 
singularities in the phase profile. Each vortex induces a 
circular velocity field about it, with the magnitude falling 
off as the inverse of the distance from the vortex center, 
and the direction determined by its charge. The motion 
of each vortex is governed by the velocity field of the 
condensate, which is typically dominated by the velocity 
field induced by nearby vortices. In our example, this 
causes each of these lines to spin in pin- wheel fashion. 
In the 17.5 ms frame, this motion has induced collisions 
between pairs of vortices (from x > and x < 0), re- 
sulting in annihilations of the vortices, with the energy 
carried off in sound waves, (seen as curved density de- 
fects) which quickly propagate towards the condensate 
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edge and break up. What remains are four vortices, as 
seen at 39 ms. At 53 ms these vortices form two vor- 
tex pairs, that, due to the small distance between each 
pair, move very quickly towards the condensate center, 
on a collision course. Upon close approach with the other 
pair, the vortices make very sudden 90 degree terms and 
switch partners. Associated with this collision is a rather 
noticeable sound pulse resulting in a ring like shape, seen 
at 57.5 ms. However, the four vortices survive and due to 
their interaction with the BEC edge, are sent around the 
edge of the BEC (79 ms) and eventually back along the 
same collision course (116 ms). This time the pairs are 
so close together that their density profiles greatly over- 
lap (compared with the first approach, 53 ms). Upon 
the second collision, at 122 ms, rather than make the 
90 degree turn, the vortices are destroyed (no phase sin- 
gularities are visible after the collision) and the density 
defects propagate through each other and continue along 
the z axis (123 ms). 

This sequence of vortex dynamics demonstrates a myr- 
iad of interesting features associated with the finite com- 
pressibility of the superfluid. In the limit of no com- 
pressibility, there is no sound emission and the vortex 
dynamics is completely conservative. In the absence of ef- 
fects associated with the background condensate density, 
a vortex pair collision would look much like the collision 
at 57.5 ms, but with no associated sound emission, and 
with the size of the outgoing vortex pairs d (at 57.5 ms), 
the same as the incoming size (53.0 ms). However, due to 
the finite speed of sound, the vortices can exchange en- 
ergy with phonons, dissipating some of their energy, as in 
the first collision of pairs at 57 ms, and in some cases an- 
nihilating, as in the second collision at 122 ms. In fact, 
the dissipation experienced during the first collision is 
what ultimately lowered the energy of each vortex pair, 
causing the pairs to to be even smaller as they entered 
the second collision. 

We note that the number of vortices created with the 
light roadblock and their subsequent dynamics can be be 
controlled both by the size of the defect (which is ulti- 
mately controlled by the number of photons in the input 
probe pulse ) an d anisotropy in the x — z plane. For exam- 
ple, in Ref. [lj it was shown how a more circular shaped 
BEC leads to the production of many more vortices in 
each soliton. This geometry would be more favorable for 
minimizing effects due to BEC density variations. A suf- 
ficiently elongated trap would produce only two per soli- 
ton, and may be more useful for tightly controlling the 
vortex trajectories and thus ensuring a direct collision is 
induced. 



D. Numerical study of dissipation in vortex pair 
collisions 

To study this collision induced dissipation more sys- 
tematically, we performed a series of numerical calcula- 
tions of vortex pair collisions in a homogenous BEC. This 



was accomplished by first relaxing a condensate in a trap 
consisting of a large square flat potential surrounded by 
steep trapping walls. The product of the homogenous 
normalized 2D density and the nonlinearity 
determine the characteristic length scale (the healing 
length) C — h/ (2mK2D\w° S> \ 2 ) 1 ^ 2 , speed of sound c = 
(k2d|V' (0) | 2 /w) 1/2 , and time scale t = h/ ^d]^? ■ 
The 2D CPE (1221) can then be renormalized into dimen- 
sionlcss form and we quote our results for this section in 
these units. Our homogenous trap region was chosen to 
by 78£ x 78£ and for numerical calculation was broken 
up into a grid with spacing 0.21£. 

After relaxation we then impose four vortices: posi- 
tive charge ones (clockwise velocity circulation) at (x = 
7A(,z = 14.7C) and (-7.4C-14.7C) an d negative ones 
at (-7.4C, 14. 7C) and (7A(, -14.70- In cach case we im- 
pose the e ±J< ^ pha se and a numerical estimate for the 
density profile |l5|. This configuration is then relaxed 
further. As it relaxes, each vortex pair shrinks (i.e. the 
vortices drift towards x = and towards each other), as 
bringing vortices of opposite charge together lowers the 
overall energy. We then stop this relaxation at various 
times and begin a real time evolution of the GPE. The 
circulations of the vortices are such that the pairs move 
towards each other and collide. In this way we can sys- 
tematically study how vortex pairs of various initial sizes 
din act upon collision. For both relaxation and propaga- 
tion of the GPE we used a Crank-Nicholson algorithm, 
with time steps O.Olio. 

Several examples of collisions are shown in Fig. [7] In 
the top example, d{ n = 7.0£ and the vortices travel along 
the paths indicated by the arrows. Upon collision the 
vortices turn 90 degrees and switch partners. There is a 
barely perceptible sound wave upon the collision. This 
sound radiation is due to the sudden acceleration experi- 
enced by the collision. In the second row d m — 5.4£, and 
a noticeably larger sound wave is generated (due to the 
larger acceleration experienced by a smaller, and there- 
fore faster, pair). Even though the density profiles of the 
pairs significantly overlap, the phase singularities remain 
in the wavefunction %p\ after the collision as the vortices 
travel out. This example is analogous to the first col- 
lision in Fig. (57 ms). In the third row, d m = 4.6C- 
In this case, the sound wave is bigger still, and the out- 
going vortices have completely merged. Examining the 
phase indicates also that the phase singularities annihi- 
late shortly after the collision. These structures bear 
resemblance to the 'lumps' studied in j5^. In the bot- 
tom case, dj n = 3.3C, the vortices annihilate early on 
in the collision, before the trajectories of the vortices 
have been significantly altered, and the remaining density 
waves simply pass through each other. This is similar to 
what happened in the second collision of Fig. 0(122 ms). 
It is interesting to observe how the sound wave energy 
emission pattern is strikingly different between the third 
and fourth cases of Fig. \7\ with the the energy emission 
primarily horizontal in the former and vertical in the lat- 
ter. We also observed an intermediate case (d[ n = 4.2£), 
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FIG. 7: Examples of vortex pair collisions, with each row 
showing two snapshots of a separate simulation with different 
initial pair distances di n . Note that here white indicates high 
density (in contrast to Fig. The left (right) hand column 
shows the vortices before (after) the collision of the pair. The 
initial distances of the pairs are, from top to bottom, dm/C = 
7.0, 5.4, 4.6, and 3.3. The arrows in the top row indicate the 
vortex direction. 



in which the sound energy propagated away with near 
perfect cylindrical symmetry. 

In Fig. Ufa) we plot a series of trajectories followed by 
the vortices during these pair collisions (due to the sym- 
metry, we only plot the vortex in the x > 0, z > quad- 
rant in each case). Initially all are following a trajectory 
almost purely in the — z direction, with velocities propor- 
tional to the arrow lengths in the figure. The input size 
of the vortex pairs c?i n is twice the initial x position. For 
sufficiently large vortex pairs d ^S> C the velocity varies as 
1 /d, in agreement with the theory of an incompressible 
fluid, however, as d becomes comparable to £ the veloc- 
ity begins to saturate, reaching a value 0.48 cq for the 
smallest pair we observed, d m = 2.18£. The phase singu- 
larities disappeared when the GP relaxed the vortex pair 
sizes to smaller than this. As each vortex approaches the 
counter-propagating pair near the origin, it begins to feel 
the pair's influence and turn to the +x direction, eventu- 
ally forming a new pair with its counterpart in the other 
pair. Again, for d ln £ vortices, the system is conser- 
vative and the new vortex pairs propagate out with the 
same size dout = ck a . However, as we saw in Fig.0 large 




2 4 6 8 10 

FIG. 8: (a) The trajectories (specifically of the phase singu- 
larities) followed by the upper-right vortex for several differ- 
ent input sizes d[ n . The arrow lengths are proportional to the 
observed velocity of the pairs at z = 10( (with the longest ar- 
row corresponding to 0.48 Co). In four of the cases the vortices 
annihilated at the locations indicated with the dots. In the 
other cases, the phase singularities survived the collision, (b) 
For a series of cases, we plot the incoming pair size d[ n (calcu- 
lated when they are at z — ±10£) versus the outgoing size di n 
(when they reach x = ±10£). (c) The energy of the vortex 
pairs versus the input distance d (black dots), calculated by 
taking half the energy of the two vortex pair configuration 
above the ground state energy. The horizontal line indicates 
the point at which the phase singularities disappear. The 
blue dots then indicate the pair energies, calculated via the 
observed d out . 



sound waves emit energy for smaller vortex pairs, and as 
a result, the outgoing pairs have lower energy and thus 
dout < d ln . The outgoing versus incoming sizes are plot- 
ted for a series of cases in Fig.|Hfb) and the trend towards 
more dissipation for smaller d- ln is clearly seen. For the 
four smallest d- m plotted in Fig. Eta), the dissipation is 
sufficient to induce an annihilation event at the locations 
indicated by the dots (thus these cases are not plotted in 

Fig. EM). 

Since the fundamental mechanism here is energy dissi- 
pation via sound emission, it is useful to analyze the nu- 
merical simulations from an energy point of view. To do 
this, we calculate the energy of our relaxed wave function 
at the beginning of each simulation. By subtracting from 
this the ground state energy (without vortices) and halv- 
ing the result (since there are two vortex pairs), we obtain 
a vortex pair energy as a function of size E pa i r (d). This 
quantity is plotted for each d[ n as black dots in Fig. |S{c) . 
The horizontal line indicates the point at which no phase 
singularity was observed (d < 2.18C). We then plot, us- 
ing the observed outgoing size d out and our numerically 
calculated energy dependence E pa i I (d), the outgoing pah- 
energy and plot it for versus d in as blue dots. The differ- 
ence between the two black and blue dots thus indicates 
the energy carried off by phonons in each case. For the 
four smallest incoming pairs, the energy dissipation was 
large enough to that the pair energy fell below the hor- 
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izontal line and the vortex phase singularity was elimi- 
nated. 



V. CONCLUSIONS 

We investigated methods for imprinting desired pat- 
terns of phase singularities on atomic matter waves by 
means of em fields inducing transitions between different 
hyperfine levels of the atom as well as by means of con- 
trolling the ultra-slow light propagation inside the atomic 
cloud. By appropriate field superpositions one may con- 
struct phase singularities around the node points of the 
em field amplitude which can then be imprinted on the 
BECs using coherently driven atomic transitions. In par- 
ticular, we showed that a simple superposition of two 
Gaussian laser beams, which exhibit a different beam 
waist, is sufficient to imprint a phase singularity forming 
a closed circular loop with a 2w phase winding around 
the loop. This provides a controlled method of creat- 
ing vortex rings on atomic superfluids, where the size, 
the position, and the orientation of the ring can be en- 
gineered by modifying the parameters of the beam su- 
perposition. Moreover, we showed that the imprinting 
techniques can be generalized for imprinting vortex ring 
pairs and multi-vortex systems which are especially in- 
teresting in the studies of vortex reconnections. The ex- 
perimental technology of diffractivc optical components, 
such as optical holograms, may be particularly useful in 
generating the light beams with the desired phase and the 
amplitude profiles. We investigated in detail one specific 
example of vortex reconnections resulting from the im- 
printed phase profile by considering two merging vortex 
rings in a trapped, two-component atomic BEC. We used 



the parameters of 87 Rb experiments and calculated the 
energy distribution in the atom cloud during the recon- 
nection process. 

In addition to direct imprinting of phase singularities 
on the BECs, we also explored the recently developed 
experimental technique of preparing defects by means of 
abruptly distorting the ultra-slow light propagation in 
atomic BECs with 'light roadblocks'. We showed that 
in a tightly-confined 2D configuration this technique is 
particularly useful in the studies of point vortex colli- 
sions. Similarly to the reconnection dynamics of the 
vortex rings in an isotropic trap, we explored the role 
of sound emission in the collision processes of the point 
vortices and showed how the binding energy of the pair 
can be dissipated during close range collisions. Suffi- 
ciently large energy dissipation leads to the interesting 
phenomenon of vortex annihilation. 

Our results point the way towards a number of exper- 
iments that can be implemented with present day cold- 
atom and laser technology, and that could allow novel 
and controlled experimental investigation of complicated 
vortex interactions. 
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